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Abstract 

Using a data sample of 282 fb _1 collected by the Belle experiment at the KEKB e + e~ collider, 
we study D° — > ir~£ + v and D° — > K~i + v decays (I = fi,e) in e + e~ annihilation. We identify 
D*+ -» D°ir + decays by using the mass of the system recoiling against a fully reconstructed tag- 
side D* or D meson, allowing for additional primary mesons from fragmentation. Using a novel 
global reconstruction method that provides very good resolution in neutrino momentum and in 
Q 2 = (Pe + Pu) 2 \ we reconstructed D° — > tt~£ + v and D° — > K~£ + v decays. From these events we 
measured the branching fraction ratios BR(D° -> irev)/BR(D° -> Kev) = 0.0809±0.0080±0.0032 
and BR(D° -> ttuu) /BR(D° -► Kfiv) = 0.0677 ±0.0078 ±0.0047, and the semileptonic form factor 
ratio U(D° -► tt'1 + v) 2 / f+(D° -► K~l + v) 2 • \V cd \ 2 /\V cs \ 2 \ q 2 =0 = 0.041 ± 0.003 ± 0.004 , where the 
errors are statistical and systematic, respectively. 

PACS numbers: 13.20.Fc,14.40.Lb,13.66.Bc 
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INTRODUCTION 



Semileptonic decays of heavy-to-light mesons are of special interest since they are well 
suited to determine CKM matrix elements such as V ub , V cd , V cs . However, some of these 
parameters — in particular V u t — are not yet measured with satisfactory precision. The 
hadronic current in the semileptonic decays of B and D mesons is parametrized through 
form factors that depend on the invariant mass q 2 of the exchanged W boson. In the 
past, the form factors for the two decays, /s(g 2 ) and fn(q 2 ), have been calculated in the 
quenched approximation of lattice QCD[1]. Results from unquenched calculations have only 
recently become available [2-4]. Imprecise knowledge of form factors is the main source 
of uncertainty in the extraction of the CKM matrix elements in semileptonic decays [5]. 
When \V u b\ is measured via B° — > 7i~£ + v, one uses the form factor /b(<7 2 ); similarly, for \V c d\ 
via D° — > ir~l + is, the form factor /z)(g 2 ) is used. In the ratio of these form factors some 
of the theoretical uncertainties cancel [6]. Accurate information on the form factor /d(? 2 ) 
from the decay D° — > 7r~£ + v can be related to the measurement of /b^ 2 ) an d hence also 
\V u b\ and can test the validity of an effective pole ansatz [7] or the ISGW2 model [8]. The 
semileptonic form factor has also been investigated by CLEO [9], FOCUS [10], BES [11], 
and CLEO-c [12]. 

In addition to the semileptonic decay D° — > tt~£ + i> we also investigate the channel D° — > 
K~£ + u, which has much less background and provides about 10 times higher statistics. From 
a study of the two semileptonic channels we can evaluate the relative branching fractions of 
these decay modes. 

DATA AND MONTE CARLO SETS 

The data collected by the Belle detector at the center of mass (CM) energy of 10.58 
GeV (T(4S f )) and 60 MeV below that, corresponding to a total integrated luminosity of 282 
fb _1 , was used. The Belle detector is a large-solid-angle magnetic spectrometer that consists 
of a multi-layer silicon vertex detector (SVD), a 50-layer central drift chamber (CDC), an 
array of aerogel threshold Cerenkov counters (ACC), a barrel-like arrangement of time-of- 
flight scintillation counters (TOF), and an electromagnetic calorimeter (ECL) comprised 
of CsI(Tl) crystals located inside a superconducting solenoid coil that provides a 1.5 T 
magnetic field. An iron flux-return located outside of the coil is instrumented to detect K\ 
mesons and to identify muons (KLM). The detector is described in detail elsewhere [13]. 
Two different inner detector configurations were used. For the first sample corresponding 
to roughly half the statistics, a 2.0 cm radius beampipe and a 3-layer silicon vertex detector 
were used; for the second sample, a 1.5 cm radius beampipe, a 4-layer silicon detector and 
a small-cell inner drift chamber were used [14]. 

Simulated Monte Carlo (MC) events were used for checks of the reconstruction method 
and background determination. The MC includes T(4S f ) — > BB and continuum (qq, where 
q — c, s, u, d) events, generated using the QQ generator [15] and processed through a 
complete GEANT-based simulation [16] of the Belle detector. We refer to this as the generic 
MC sample. We also use MC samples of signal events for optimisation of selection criteria 
and signal reconstruction efficiency determination. 
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ANALYSIS STRATEGY AND EVENT SELECTION 



To measure semileptonic decays D° — > h~£ + u (h = n or K, £ = e or //), we consider 
event topologies of the type e + e _ — > D^D*X, where the system is referred to as 
the tag side, the D* system is referred to as the signal side, and X denotes additional 
mesons including 7r°, 7r =l= and K ± (Fig. 1). We use a global reconstruction technique that 
requires all the final state particles in the event, with the exception of the neutrino from 
the signal D semileptonic decay, to be detected. In this paper, the reconstruction of charge 
conjugate modes is implied throughout. We found in agreement with [17] (which studied 
the case of no additional particles) that the e + e~ — > DDX cross section is negligible, while 
for DD* and D*D* the cross sections are so large that in spite of decreasing reconstruction 
efficiencies with increasing number of fragmentation particles, statistics can be significantly 
enhanced by our reconstruction method. Instead of reconstructing decay chains step-by- 
step by applying intermediate selection criteria, all possible decay chains matching certain 
patterns (event topology and loose mass windows) are considered, leading to a possibly 
large number of intermediate combinations in a single event that is pruned at later steps in 
the reconstruction chain. This strategy in many cases automatically solves the problem of 
ambiguities that would arise in a sequential reconstruction, and also allows looser selection 
criteria for individual particles, thus giving a higher reconstruction efficiency Appropriate 
selection criteria, described below, are applied to select D° semileptonic decays. The case of 
multiple signal candidates is treated by assigning each remaining candidate the same weight. 

Tracks are detected with the CDC and SVD. They are required to have at least one 
associated hit in the SVD and an impact parameter with respect to the interaction point in 
the radial direction of less than 2 cm and in the beam direction of less than 4 cm. Tracks 
are also required to have momentum in the laboratory frame greater than 100 MeV/c. A 
likelihood ratio for a given track to be a kaon or pion was obtained by utilising specific 
ionisation energy loss measurements made with the CDC, light yield measurements from 
the ACC, and time of flight information from the TOF. Lepton candidates were required to 
have momentum larger than 500 MeV/c. For electron identification we use position, cluster 
energy, shower shape in the ECL, combined with track momentum and dE/dx measurements 
in the CDC and hits in the ACC. For muon identification, we extrapolate the CDC track 
to the KLM and compare the measured range and transverse deviation in the KLM with 
the expected values. Photons are required to have energies in the laboratory frame greater 
than 50 MeV. Neutral pion candidates were reconstructed using 77 pairs with invariant 
mass within ±10 MeV/c 2 of the nominal 7r° mass. These candidates are also subject to a 
mass constrained vertex fit, assuming they arise from the interaction point, for which the 
confidence level is required to be greater than 0.1. 

We reconstruct tag-side D ± and D° candidates in the channels — > K ± mr and D° — > 
K ± rni where the pions can be charged or neutral, and n — 1, 2, 3. Depending on the channel 
the candidate mass window varies in the range 5 — 60 MeV/c 2 around the nominal D mass, 
corresponding to ±5(7. These candidates undergo a mass constrained vertex fit and we 
retain only those for which the fit is successful. An attempt is then made to reconstruct D* 
candidates on the tag side via D* + -> D°n + , D* + -> D + ir°, D*° -> D°n° , D*° -> and 
charge conjugate modes. The mass window around the nominal D* mass is 5 (20) MeV/c 2 
for the non-radiative (radiative) mode. We require a successful mass constrained vertex fit 
of the D* candidate. At this stage of the reconstruction we either have a fully reconstructed 
jj±o* or jj±o candidate on the tag side. 
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FIG. 1: Scheme of reconstructed event topology. Left: DD* events- D(tag)D* ± (signal). Right: 
D*D* events-D*(tag)£)* :t (signal). 

Along with possible additional primary mesons, we search for signal side D* ± in the recoil 
of the tag side: the 4-momentum of the candidate G»o*±cand) is reconstructed using the 
4-momenta of the beams (pbeam), that of the tag-side Pu(*) and of the primary meson system 
Px as p.o*±cand = Pbcam — Pd(*) ~ Px- The accuracy of the D* ± candidate momentum is 
improved by applying a mass constrained vertex fit to the recoil particle. For candidates 
where the fit to D* has a confidence level probability larger than 0.1%, an analogous fit is 
repeated for a recoil D° from D*^ — > D°7rf decay. The charge of the slow pion tt s is required 
to be the same as that of the kaon arising from the tag side D. The small available phase 
space for D* + — > D°n + in addition to the mass constraint leads to a very sharp peak in the 
D° candidate mass distribution (Fig. 2). Finally, neutrino candidates in the semileptonic 
D° decay are identified on the basis of the missing mass squared ^ miss , where the missing 
4-momentum is calculated as 

Pmiss P_D,signal— side PK/tt Pe/fji- (1) 

Neutrino candidates are retained if their CM energy is greater than 100 MeV, if the remaining 
photon candidates not used in the reconstruction have an energy sum in the CM frame of 
less than 700 MeV and finally if the charges of the lepton and the slow pion are the same. 

The signal region for D° semileptonic decays is defined as |m m i SS | 2 < 0.05GeV 2 /c 4 . Our 
method yields a good resolution on the neutrino mass and the momentum transfer q 2 by 
applying a mass constrained vertex fit for neutrinos. The latter is found to be a q 2 = 
0.0145 ± 0.0007 GeV 2 /c 2 in MC signal events. Although in the fit particles by default 
are assumed to originate from the interaction region, the influence of the D° decay distance 
on the reconstructed q 2 value has been corrected for by shifting the D° vertex in the direction 
of the fitted D° momentum by a distance corresponding to the mean lifetime of the meson. 
No significant change of the q 2 value or resolution has been observed. 

For future reference the multiplicity of an event is defined as the total number of 7r ± , tt° 
and K ± mesons not assigned as decay products of the tag-side D or signal-side D*. 
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After the selection 75% of remaining events have one candidate while a negligible fraction 
have greater than three candidates. Multiple candidates in events are mainly due to an 
interchange of particles within a particular decay chain, each is assigned equal weight within 
an event. A possible bias in the measurement of q 2 that may arise due to events where the 
lepton and meson are interchanged, a double mis-assignment, was checked with candidate 
D° — > K~l + v events, and found to be negligible. 

BACKGROUND STUDIES AND BACKGROUND SUBTRACTION 
Non D°-events and events with badly reconstructed D° 

Fig. 2 (top) shows the composition of the inclusive signal side D° invariant mass spectrum 
for generic MC, where the reconstruction of the neutrino candidate has yet to be performed. 
The generic MC sample consists of events that may contain D°, namely T(4S) — > BB 
(bottom), e + e~ — > cc events (charm), and a pure background sample of e + e~ — > qq where q = 
u, d, s (uds). Signal events have been distinguished according to MC generator information. 
Most of the background comes either from uds or from a misreconstruction of the D°. There 
is also a small contribution from the i?-meson samples (~ 1.4%). Performing the full signal 
side reconstruction and applying all signal selecting requirements as described above reduced 
the relative background level to approximately 1% (10%) for D° — > K~l + v (D° — > 7r~£ + z/), 
with the remaining background dominated by charm events with an incorrectly reconstructed 
D°. Several constrained fits that are performed in the reconstruction may make background 
events look more like signal (peaking backgrounds). This possibility has been investigated 
in the generic MC sample and no peaking of the background has been observed. 

The background is estimated by taking the D° mass spectrum from a wrong-sign D° 
sample defined by opposite charges of the kaon from the tag-side D decay and the signal- 
side slow pion. The right- and wrong-sign D° mass spectra in generic MC are shown in Fig.2. 
Subtracting the wrong-sign D° mass distribution from the right-sign distribution results in 
almost complete cancellation of the background in the latter, implying that the non-D° 
background is dominantly charge-uncorrelated [23]. The bottom plot in Fig. 2 shows that 
D° candidates are also found in the wrong sign sample, with a yield approximately 10% of 
that in the right-sign sample. 

Since the charge-correlated non-_D° background is very small, only the charge- 
uncorrelated part is considered, and a systematic error is assigned for the charge-correlated 
part, defined as the difference of background in right- and wrong-sign, equal to its relative 
contribution as observed in MC. To subtract the small component of real D 0, s in the wrong- 
sign sample, we use the difference of the right- and wrong-sign sample as model of the true 
D° signal shape, and fit this to the wrong-sign sample. For this purpose, the wrong-sign 
background is approximated as a second order polynomial in the signal region, 1.862 — 1.867 
GeV/c 2 . The fitted wrong-sign D° component is subtracted, and the remaining distribution 
is used to represent the shape of the non-.D° background. It is normalized to the data in the 
sideband 1.84 — 1.85 GeV/c 2 of the D° mass spectrum thus obtaining the amount of non-£> 
background in the signal region (which is defined as the region selected by the criterion on 
the confidence level of the D° mass-constrained vertex fit). This method is confirmed with 
generic MC (using D° — > 7r~£ + z/ selection criteria): the measured background is 6.5 ±1.4 
events, compared to the true value of 6. 

Yields of the measured non-_D° background are given in Table I. The systematic error 
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FIG. 2: Inclusive D° mass spectrum in generic MC events for right- (top) and wrong-sign (bottom) 
samples. The invariant mass shown is that prior to the D° mass-constrained fit. 

for this measurement consists of a statistical part (from background shape and sideband 
statistics), a part due to the subtraction of wrong-sign signal events, the effect of the small 
component of charge-correlated background measured with MC, and of possible bias of 
the background shape due to selection criteria, which has been estimated by comparing 
the measurement with the result using loose selection criteria, i.e., only those selecting D° 
(without the requirement on the confidence level of the mass constrained fit). Individual 
uncertainties are listed in Table II. 



Background from semileptonic decays with a misidentified or undetected meson 

The background from real semileptonic decays, with an incorrectly identified meson or 
additional mesons lost in reconstruction, is highly suppressed by the good neutrino mass reso- 
lution (see measurements below). For the signal channel D° — > 7r~£ + z/, there are backgrounds 
from D° -> K~£ + u, D° -> p-£ + u (with p~ -> t^tt ), D° -> K*-i + u (with K*~ -> K°ir-), 
as well as non-resonant decays into the latter two final states. The most significant back- 
ground is due to D° — > K~ t+v since the branching fraction is approximately 10 times larger 
than the signal. 
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For the signal channel D° — > K~£ + v, it was verified using MC that backgrounds from 
D° — > 7r~£ + z/ and D° — > p~£ + v are completely negligible due to the much smaller branching 
fractions. The latter channel is suppressed even further by the selection criteria, due to the 
additional missing momentum taken by the ir° from p~ — > The channel D° — > K*~£ + v 

contributes a small amount via — > K~tt°. 

To measure the semileptonic background in the -D — > 7r~£ + z/ sample from D° — > K~£ + v, 
the distribution of simulated -D — * K~£ + v decays that pass the selection criteria was nor- 
malized to the signal peak of the D° — > K~£ + v data events, then multiplied by the expected 
misidentification probability for kaons to be identified as pions. The misidentification rates 
were corrected for known differences with real data, measured with a dedicated D* sample 
in bins of meson momentum and polar angle. The remaining contribution of D° — > K*~£ + v 
and D° — > p~£ + v was obtained from MC and normalized in the region m 2 u > 0.3 GeV 2 /c 4 . 
For D° — > K~£ + u, the only semileptonic background stems from D° — > K~*£ + u, which was 
determined with normalized MC as described above. 

The systematic error in the estimate of the semileptonic background includes uncertainties 
due to MC statistics, in the misidentification probability correction and in the ratio of K* 
and p contributions. Yields of the semileptonic background are shown in Table I; their 
contributions to the systematic error are listed in Table II. 

Hadronic D°-decays with fake leptons 

As the fake rate for muons is about an order of magnitude larger than that for electrons, 
the background from hadronic -D°-decays with fake muons is much more significant. To 
study and further categorize this kind of background, it is useful to look once more at 
wrong-sign data, in this case the lepton has opposite charge to that of the slow pion. The 
large D° — > K~n + channel can contribute to the background for D° — > K~£ + v (right-sign), 
and D° — > tt~£ + v (wrong-sign) [24]. The wrong mass hypothesis allows some energy for 
a possible neutrino candidate, but with the additional requirement E u > 100 MeV this 
background is removed for both cases. However, in the decay D° — > K~it + it the it may 
have enough momentum to fake a neutrino that passes the E v selection criterion. For 
D° — > K~£ + v ( D° — > 7r~£ + z/), this background contributes to the right-sign (wrong-sign) 
sample. The background channel D° — > 7r~7r + 7r° contributes equally to right and wrong sign 
samples. 

To measure the shapes of these backgrounds, special samples are used wherein identified 
kaons and pions are intentionally misidentified as leptons. As these raw shapes also contain 
some amount of non-D° background, the same procedure as described in the previous section 
is used to measure and subtract this part. A linear combination of the shapes of intentionally 
misidentified kaons and pions is fitted to the measured m 2 u distribution of the wrong-sign 
sample (see Fig. 3). Assuming that the K ± and ti^ misidentification rates do not depend 
on the charge of the slow pion, the same amount of background contributes to the wrong- 
sign and the right-sign samples. Since the statistics for wrong-sign data samples in the 
D° — > K~£ + v channel is very small, we use the measured fake-rates in D° — > ir~£ + v for the 
kaon modes as well. Numerical results can be found in Table I. 

It has been checked using the generic MC sample that the shape of the background is 
relatively insensitive to the differences in the kinematics of correctly identified and mis- 
identified mesons; the possible small differences have been included in the systematic error 
(see Table II). The statistical error in the fits to the wrong-sign sample is included in the 
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D° -> 7iev 



D° -> 7IUV 




FIG. 3: Fitted hadronic background shapes for the pion channels in the wrong-sign data sample 
(defined by opposite charges of the slow pion and the lepton); the two prominent peaks correspond 
to background from D° — > K~ir + 7r° (lower m^ liss ) and D° — > K°it + it~ (higher m^ iss ) 

systematic uncertainties as well. 

Fig. 4 shows the measured backgrounds and the signal for all four semileptonic decays 
of the D°. Table I gives the yields for the background sources. The MC sample was used 
to check that all the significant background components are taken into account using the 
above described methods. 

BRANCHING RATIOS OF SEMILEPTONIC DECAYS OF D° 

After background subtraction, the yields for the four studied semileptonic D° decay 
modes are given in Table I. The contributions to systematic uncertainties for each individual 
background source are listed in Table II. 

The dependence of signal yields on the multiplicity has been studied with generic MC. 
Although the efficiency decreases as the multiplicity of the event increases, the ratio e v /e K , 
where e n (e K ) is the efficiency of the D° — > (K~£ + u) channel, is flat (Fig. 5, jndf = 

0.84). The observed multiplicity in the data is slightly different from the simulated one. In 
order to estimate the systematic error due to this difference, we reweighted the simulated 
efficiencies according to the data multiplicity distribution. Since the D° — > ty^£ + u and 
D° — > K~£ + v decay samples are topologically very similar, we used a sample of events 
where only the selection on the signal-side D° meson was applied. The efficiency ratio 
for the reweighted inclusive D° sample and the D° — > K~l + v sample compared to the 
simulated ratio yields a correction factor of 0.981 ± 0.028. Since this is compatible with 
unity, no correction is applied to the efficiency ratio; instead, the uncertainty of this factor 
is included in the systematic error. The values for the efficiencies, corrected and averaged 
over multiplicity, are listed in Table III. 
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FIG. 4: Decomposition of background for the four semileptonic channels. In each case the lower 
plot shows the background-subtracted signal. 

Using these efficiencies, the relative branching ratios are obtained and are shown in Table 
IV. We have divided our data sample into low- and high-multiplicities sets of three or fewer 
and four or more pions (kaons), respectively, and found consistent results. The results are 
also in agreement with recent measurements from CLEO [9, 12] and FOCUS [20] 
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DISTRIBUTION 



The square of the four-momentum transfer q in the semileptonic channel is given by 
I 2 — (Pu + Pe) 2 , where p is the four momentum of the specified particle. The good neutrino 
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channel 


K~e+v 


K~n+v 


7T e + u 


TT fl + V 


total yield 


1349 


1333 


152 


141 


signal yield 


1318 


1248 


126 


107 


statistical error 


37 


37 


12 


12 


total systematic error 


7 


24 


2.6 


6 


non-D° background 


12.6 ±2.2 


12.4 ±4.8 


12.3 ±2.2 


11.5 ±4.5 


semileptonic background 


6.7 ± 2.6 


10.0 ±2.5 


11.7 ± 1.2 


12.6 ± 1.9 


hadronic background 


11.9 ±5.6 


62.1 ±23.9 


1.8 ±0.7 


9.7 ±3.7 



TABLE I: Final signal yields and backgrounds. 



channel 


K-e + v 




7T e + v 




non-D° background: 




statistics of background shape 


4.2% 


8.1% 


6.6% 


13.4% 


subtraction of wrong sign 


1.0% 


1.0% 


1.3% 


1.3% 


charge correlated background 


6.3% 


13.1% 


5.5% 


11.9% 


shape minimum bias 


15.7% 


35.2% 


15.3% 


34.1% 


sum in quadrature 


17.5% 


38.4% 


17.6% 


38.5% 


semileptonic background: 




statistics MC samples 


39% 


25% 


9.7% 


14.4% 


fake rates ratio MC/data 


4% 


3% 


2.9% 


3.2% 


uncertainty of K~*t r v / p~ l^v 


< 1% 


< 1% 


1.4% 


1.7% 


sum in quadrature 


39.2% 


25.2% 


10.2% 


14.8% 


hadronic background: 




fit of wrong sign shapes, mis-id kaons 


2% 


1% 


3% 


3% 


fit of wrong sign shapes, mis-id pions 


33% 


16% 


35% 


23% 


bias of background shapes 


33% 


35% 


11% 


30% 


sum in quadrature 


46.7% 


38.5% 


36.8% 


37.9% 



TABLE II: Relative systematic errors on the amount of respective backgrounds. 

resolution leads also to a very good resolution for q 2 , described by a double Gaussian with a 
common mean of 0.00061 ± 0.00031 GeV 2 /c 2 and widths of a x = 0.0109 ± 0.0005 GeV 2 /c 2 , 
cr 2 = 0.0401 ± 0.0026 GeV 2 /c 2 , the fraction of the wider Gaussian being 0.14 ± 0.02. No 
unfolding of the q 2 distributions is needed since these resolutions are much less than the bin 
widths used. The dependence of q 2 resolution on different multiplicities is small; a(q 2 ) lies in 
the range 0.0095 to 0.0121 GeV 2 /c 2 (with a ~ 10% uncertainty) for pion/kaon multiplicities 
between and 8. The background subtraction in the determination of the q 2 distribution 
is performed using the method described earlier. The results are shown in Fig. 6 and are 
consistent with MC. 
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channel 


K~e + v 




efficiency 


(0.369 ±0.005 ±0.010)% 


(0.350 ±0.005 ±0.010)% 


channel 






efficiency 


(0.436 ±0.006 ±0.012)% 


(0.443 ±0.006 ±0.012)% 



TABLE III: Efficiencies for semileptonic decay channels, averaged over multiplicity; errors shown 
are due to MC statistics and multiplicity dependence, as described in the text. 
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0.5 
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Multiplicity 



FIG. 5: Efficiency ratio for D° — ► it £ + u/D° — ► K £ + v generic MC events as a function of 
multiplicity (number of additional mesons). The solid line is the result of a fit to a constant. 



relative BRs 



this analysis 



PDG [5] 



theor. pred. [22] 



BR(D"^K~e+u) 











e+ 




BR(D°^7T- 




•0 


BR(D°^7T- 




") 








BR(D»^7r- 


>H 








*■») 



1.002 ±0.041 stat ±0.048 syst 
1.20±0.18 stat ±0.09 syst 
0.0809 ± 0.0080 stat ± 0.0032 syst 
0.0677 ± 0.0078 stat ± 0.0047 syst 



1.12 ±0.08 



0.101 ±0.018 



1.03 



1.02 



0.086 



0.087 



TABLE IV: Relative branching ratios 
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FIG. 6: q 2 distributions for D° -» R-e+v (upper left), Z?° -► K~ ' (upper right), D° — > 7r e + z/ 
(lower left) and D° — > 7r~/i + z/ (lower right), with contributions of individual background sources 
shown as histograms. 

SEMILEPTONIC FORM FACTOR FOR D° -> tt-^+i/ AND D° -> K~£+u 

In general, for P — > P transitions (P denoting a pseudo-scalar meson), two form factors 
f+(q 2 ) and f~(q 2 ) are needed [7]. In the case of semileptonic decays, the matrix element 
is dominated by the form factor / + ; the influence of /_ is suppressed kinematically to a 
negligible (small) level for electrons (muons) by the smallness of mf compared to the scales 
of m 2 D , p 2 and q 2 . Up to order m 2 one has 



dT(D° -> G 2 F \V CC 



|2 



3? " ^M^ti < 2 > 

and an analogous expression for D° — > K~£ + i>, where p n is the magnitude of the pion 3- 
momentum in the D° rest frame [25], determining the kinematical weight of the form factor. 
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These form factors have been calculated recently in unquenched lattice QCD [2], [4]. Earlier 
calculations were done in quenched lattice QCD [1], which works best in the region of high 
q 2 where measurements unfortunately suffer from poor statistics. In the pole model [7], the 
form factor / + is described as 



pole y 



W) = n 7/ 2 s ( 3 ) 

(1 - q 2 /m 2 ) 



with the pole masses m D * = 2.11 GeV/c 2 (for D° -> K~l + v) and m D * = 2.01 GeV/c 2 (for 
Within the modified pole model [22], the form factor is given by 

The ISGW2-model [8] predicts the following expression for the form factor: 

f+(<LJ /+(o)(i + 
/+w J " (i - M? 2 - <?LJ) 2 " (i - «/(<? 2 - <?LJ) 2 1 J 

where (^ aa . is the kinematical limit of g 2 and ai is a parameter of the model. 

The normalized, bin-by-bin background subtracted and efficiency corrected q 2 distribu- 
tions are shown in Fig. 7 for different decay modes; error bars indicate the sum in quadrature 
of statistical and systematic uncertainties. We fit the measured distributions to the predicted 

i t/r 

normalized differential decay width — — of different models described above. Binning ef- 

T dq 2 

fects are accounted for by averaging the model functions over individual q 2 bins. Parameters 
/(0) and fiq^ax) cancel out in the normalized decay width. 

The fit results for the simple pole model with the pole mass as free parameter are shown 
in Fig. 7, left. The fitted pole masses are 

m P oie(D° -> K~e + v) = 1.88 ± 0.06 ± 0.03 GeV/c 2 (6) 

m po i e (D° -> K~fi + u) = 1.77 ± 0.04 ± 0.03 GeV/c 2 (7) 

m pole (D° -> n-e + v) = 2.01 ± 0.13 ± 0.04 GeV/c 2 (8) 

m pole (D° -> tt> + z/) = 1.92 ± 0.09 ± 0.04 GeV/c 2 (9) 

where the first error is statistical and second systematic (discussed below). The xV n df 
values for fits using different models are shown in Table V. Values of m vo i e {Klv) are several 
standard deviations below the mn*. The pole mass for the i^tv decay agrees within errors 
with the predicted value, m D *. The fitted pole masses are also in agreement with results 
from CLEO [9] and FOCUS [10]. 

A fit to the modified pole model is also performed, where m p is fixed to the theoretical 
pole rriD* . The results for the fit parameter a„ are: 

(s) 



a p (D° -> K~e + u) = 0.40 ± 0.12 ± 0.09 (10) 

a p (D° -> K~n + v) = 0.66 ±0.11 ±0.09 (11) 

a p (D° -> Ti-e + v) = 0.03 ±0.27 ±0.13 (12) 

a p (D° -> 7T>+z/) = 0.19 ±0.32 ±0.16 (13) 
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D° -> Klv D° -> Klv 




q 2 (GeV 2 /c 2 ) q 2 (GeV 2 /c 2 ) 



FIG. 7: Fraction of events in different g 2 -bins for D° — > K i + v and D° — > tt l + u compared with 
the fits of the simple pole model (left) and ISGW2 (right) 

Significantly non-zero values of a v {Kiv) suggest further contributions to the form factor, 
apart from the single pole at m^. 

Finally, a fit of the parameter on in the ISGW2 model yields (Fig. 7, right) 



a/(D° —> K~e + v) = 0.37 ± 0.03 ± 0.02 GeV" 2 c 2 (14) 

ai (D° -> K~fi + u) = 0.44 ± 0.03 ± 0.02 GeV" 2 c 2 (15) 

a I (D° -> n-e + u) = 0.35 ± 0.08 ± 0.04 GeV" 2 c 2 (16) 

a T {D° -> = 0.40 ± 0.10 ± 0.03 GeV" 2 c 2 . (17) 



The x 2 °f the fit to the ISGW2 model is comparable to that for the two other models (c.f. 
Table V). The prediction of [8] for the Klv mode is aj = 0.47, which is within one standard 
deviation of the K[iv result, and almost three standard deviations from the Kev result. 

To study systematic uncertainties in different fits, each fit was repeated 50000 times, 
with the individual background normalizations varied within their Gaussian errors in a 
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X 2 /ndf 


K~e+v 




7T e + v 


TT 


simple pole model 


1.08 


1.40 


0.37 


0.85 


modified pole model 


1.05 


1.35 


0.37 


0.89 


ISGW2 


1.02 


1.33 


0.35 


1.09 



TABLE V: x 2 per degree of freedom for fits to various models 



correlated manner for all q 2 bins. The mean differences from the results of default fits were 
taken as systematic errors. A similar procedure was adopted to estimate the error due to 
the uncertainty on the q 2 dependent efficiency The systematic error is included in the fit 
results above. 

The absolute (i.e. not normalized) partial decay width is related to observable quantities 

by 

dT(D° -> (ir/K)-£ + v) e K BR K N{q 2 ) 



dq 1 



r D e a 2N K 



(18) 



where ex , BRk and Nk are the average efficiency, branching fraction, and number of 
reconstructed signal events in the D° — > K~£ + u mode, respectively. The world average 
value [5] is used for the branching fraction, as well as for the D° lifetime t d . e q 2 and N(q 2 ) 
are the efficiencies and signal yields in individual q 2 bins. A comparison of the measured 
values with lattice calculations [2] is shown in Fig. 8 (left). The uncertainty of the lattice 
results is mainly systematic (errors given in [2] are statistical only; relative systematic errors 
of about 10% were added in quadrature according to [3]). The lattice results in Ref. [4] are 
very similar to those in Ref. [2], and are therefore not shown here. 

From the absolute partial decay width, also the absolute form factor can be extracted 
using Eqn. 2. Comparison of the measured values and lattice calculations [2] is shown in 
Fig. 8 (right). In order to make a more direct comparison, f+(q 2 ) is also calculated at the 
same q 2 values where the prediction is available. For interpolation we use the modified pole 
model (also used in [2]) with a p obtained from the fit to data described above. 

Using the integrated form of Eqn. 2 we determine the ratio of the form factors for Klv 
and ntv decay modes at q 2 = 0. For the q 2 dependence of the form factor we use the 
modified pole model with the values of parameter a p in Eqn. 10-13, with the weighted 
average of K and n modes calculated taking into account correlated systematic errors. A 
similar averaging procedure is done for the measured ratios of branching fractions (Table 
IV). We measure the ratio 



7T-£+v) 2 \V cd \ 2 



U(D°^K-i+isy\v cs \ 2 



q 2 =0 



0.041 ± 0.003 stat ± 0.004 syst 



(19) 



which is consistent within errors with the model-independent result using only the data in 
the first iriis q 2 bin (q 2 < 0.3 GeV 2 /c 2 ). A recent theoretical prediction for the ratio [2] is 
0.040 ±0.004. 

This result is in fair agreement with those from CLEO [9] and FOCUS [20], which measure 
slightly lower values. 
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D° -> Klv 



D° -> Klv 




q 2 (GeV 2 /c 2 ) q 2 (GeV 2 /c 2 ) 



FIG. 8: Left plot: comparison of the absolute partial decay width dT/dq 2 with a lattice calculation 
(circles); right plot: comparison of the measured form factor (error bars) and unquenched lattice 
results (circles). To make a more direct comparison, f+(q 2 ) is also calculated at the same q 2 values 
where the lattice prediction is available (shifted slightly higher to avoid overlap), using the modified 
pole model with a p parameter from the fit to the data (crosses) 



SUMMARY AND CONCLUSIONS 



In the present measurement we use fully reconstructed e + e~ — > D\*a g D* s f gnal nh events 
(h = 7r ±0 ,i^ ± ), where the tag-side decay channels were D* — > Dn^D^i and D — > K ± nir, 
n = 1 — 3 and the signal side was reconstructed from the recoil as D*^ — ► D ^ and 
D° — > K/n~£ + b>, I = e,/i. Mass constraints at each reconstruction step enabled a very good 
resolution for the momentum transfer q 2 of a q 2 ~ 0.01 GeV 2 /c 2 . 

We determined relative branching fractions 

BR(D° -> 7T-eV) , , 

BR(D o _> ^- e+ 4 = °- 0809 ± °- 008 °stat ± 0.0032 syst (20) 
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BR ( D ° ^ n = 0.0677 ± 0.0078 stat ± 0.0047 svst (21) 

BR(L>° -> K-fi+v) stat syst V ; 

in good agreement with other measurements. The normalized measured q 2 distribution was 
fitted to different models of form factors. In the Kiv decay mode, where the sample of 
events with higher statistics is selected, some significant deviations from the predictions of 
the simple pole and ISGW2 model are observed. The results of the fit to the pole masses are 
m polc {K-£ + u) = 1.81 ± 0.03 ± 0.02 GeV/c 2 and m polc {n-£ + u) = 1.95 ± 0.07 ± 0.03 GeV/c 2 . 
Using the measured values of the branching fraction ratio and the modified pole model with 
fitted parameter a p , we obtain 

f + (D° - n-£ + u) 2 /f + (D° - K-tuf ■ \V cd \ 2 /\V cs \ 2 \ q2=0 = 0.041 ± 0.003 stat ± 0.004 syst . 

(22) 

in agreement with recent measurements [9, 20]. The form factors f(q 2 ), obtained using the 
measured dT/dq 2 distribution and the world average of the Br(D° — > K£u), are in nice 
agreement with the predictions of the recent unquenched lattice QCD results [2]. 
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